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Abstract 



Distributions derived from non-extensive Tsallis statistics are closely con- 
nected with dynamics described by a nonlinear Fokker-Planck equation. The 
combination shows promise in describing stochastic processes with power- 
law distributions and super diffusive dynamics. We investigate intra-day price 
changes in the S&;P500 stock index within this framework by direct analysis 
and by simulation. We find that the power-law tails of the distributions, and 
the index's anomalously diffusing dynamics, are very accurately described 
by this approach. Our results show good agreement between market data, 
Fokker-Planck dynamics, and simulation. Thus the combination of the Tsal- 
lis non-extensive entropy and the nonlinear Fokker-Planck equation unites in 
a very natural way the power-law tails of the distributions and their superdif- 
fusive dynamics. 
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I. INTRODUCTION 



In anomalously diffusing systems, a mean-square displacement scales with time accord- 
ing to a power law, t", with a > 1 (superdiffusion) or a < 1 (subdiffusion). (The case 
a = 1 corresponds to normal diffusion.) Anomalous diffusion has been observed in systems 
as widely varied as plasma ffow [|l|, surface growth 0, and financial markets 0. A general 
framework for treating superdiffusive systems is provided by the nonlinear Fokker-Planck 
equation, which is associated with an underlying Ito-Langevin process [^H^. This in turn 
has a very interesting connection to the nonextensive entropy proposed by Tsallis the 
nonlinear Fokker-Planck equation is solved by time-dependent distributions which maximize 
the Tsallis entropy This unexpected connection between thermostatistics and anoma- 
lous diffusion gives an entirely new way to approach the study of dynamics of anomalously 
diffusing systems. In this paper we use this viewpoint to address the dynamics of financial 
markets. 

Several financial markets' indices as well as their member stocks are characterized by 
price changes whose variances have been shown to undergo anomalous (super) diffusion 
under time evolution [0,^11]. Moreover, the probability distributions of the price changes 



have power law tails 0,[T2|. An open long-term question is how best to describe these 
distributions and their time evolution. Earlier work has shown that the power-law tails can 



be described by a distribution which maximizes the Tsallis entropy [T^|Tj]. The connection 
mentioned above then suggests that the market dynamics might be controlled by a nonlinear 
Fokker-Planck equation. In this description the power-law tails and the anomalous diffusion 
arise together quite naturally. Here we show that this approach seems to accurately describe 
high-frequency intra-day price changes for the S&P500 index. It may also be suitable to 
describe the superdiffusion and power law behaviors observed in a broad range of markets 
and exchanges [Q. 

The description we will use was developed in the general context of anomalously diffusing 
systems by Tsallis and Bukman and Zanette |^ . Here we briefly summarize their results, 
which are based on a maximization of entropy subject to certain constraints. The central 
ingredient is a time-dependent probability distribution P{x,t) of a stochastic variable x. 
The point of departure from normal maximum entropy approaches is the feature that the 
entropy used is the non-extensive Tsallis entropy, 

Sg = — — (i- f p{x, ty dx) . (1) 



1-q 

The Tsallis parameter q characterizes the non-extensivity of the entropy. In the limit q ^ 1 
the entropy becomes the usual logarithmic expression S=-/ Pin P. 

Associated with the non-extensive entropy is the use of the constraints 



P{x,t)dx = l, (2) 
{x - x{t))g = J[x- x{t)]P{x, ty dx = 0, (3) 
{{x - x{t)y), = ![x- x{t)fP{x, ty dx = aq{tf. (4) 



The first of these is simply the normalization of the probability. However, in the latter two 
equations the probability distribution function is raised to the power q. Unless q = 1 these 
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are not the usual constraints leading to the mean and variance, and cXg (the 'g-variance') is 
not the ordinary variance. Notice that the Tsallis nonextensivity parameter q is independent 
of time. 

Maximizing the Tsallis entropy subject to these constraints, for fixed g, yields 

p(x, ^) = ^ {i + m{q - - m'}'"^' ■ (5) 

Here parameters Z (a normalization constant) and (3 are Lagrange multipliers associated 
with the first and third constraints, and are given by 

(l _L 1 

z{t) = ^^ (6) 

^^^^ ^ 2a^{tyZ{ty-^' 

where B{x, y) = T{x)T{y)/T{x + y) is Euler's Beta function. 
The ordinary variance of the distribution Eq. is 

r 1 ^< 5 



a\t) = ((x - x{t)f)^ = (5-3g)/3W ' ^ ^ I (8) 

Hence for applications to data of finite variance, the Tsallis parameter q must lie within the 
range 

1 < g < 5/3. (9) 

The Tsallis function Eq. can be viewed as a least biased probability distribution compat- 
ible with observed stochastic data with a certain mean x{t) and variance cr^(t) (but with a 
particular choice of q). 

An important property of the probability distribution function Eq. (|^) is that, with 
appropriate time-dependent parameters, it is the solution of a time evolution equation which 
leads naturally to anomalous diffusion . Consider the nonlinear Fokker-Planck equation 

where F{x) is a linear drift force, F{x) = a — bx. One can show that a function of the 
form Eq. (^) solves this as long as 

q = l + fi-u (11) 
and the time dependences of the parameters are given by 

- ^ "^^'2^^^ + ^^DP{to)Z{tof'' - bZ^+^ = 0, (12) 
Pit) _ (Zito)]'' 

P{to)-[z{t)) ' ^^'^ 
dx 



By comparing Eqs. and ([T3|) one easily shows that the nonhnear Fokker-Planck equation 
preserves the norm (/ P{x,t)dx) only if /i = 1. Hereafter we specialize to this case. Then 
Eqs. ([Tl||g|T|) give 

Z{t) is then obtained from Eq. (^). 



Pity 
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II. APPLICATION TO A MARKET INDEX 

It is apparent from the results sketched in the previous section that the Tsallis probability 
distribution function P{x,t) has properties which make it a good candidate for describing 
the anomalous diffusion of financial market indices, member stocks, and currency exchanges. 
First, in the large x limit P{x,t) becomes a power law distribution, 

P(x,t)~a;"^. (16) 
This is in keeping with the observed power law tails in market distributions and, indeed. 



certain market distributions have been well fit by a distribution of the Tsallis form |T3| , P 
In fact this fit to market data has a longer history — the Tsallis function is an extension to 
continuous values of q of the student-t distribution, which has been known for some time to 



provide a good fit to certain market data [|15 



What has not been widely recognized in the finance-related applications of statistics is 
that the Tsallis distribution solves the nonlinear Fokker-Planck equation. Consequently it 
may provide a framework for understanding the dynamics of certain market data, including 
anomalous diffusion. It is this possibility that we will now discuss. 

For application to market indices the Tsallis parameter q must lie within the range 
Eq. (y), ensuring that the regular variance remains finite and evolves for moderately long 
times as 

cT^{t) ~ l/(3{t) ~ (17) 

[using Eq. (|^JT^) in the case 6 <^ 1] . Thus the Fokker-Planck equation can describe superdif- 
fusive processes (with 1 < g < 5/3). 

Here we investigate one test case, the S&P500 stock market index, using fairly high- 
frequency 1-minute- interval data collected from July 2000 to January 2001. This data con- 
sists of a set of prices (index values) p{t) at discrete times r (r = 1,2,..., 50000 trading 
minutes) obtained from the Terra-Lycos QCharts server |jT6|. (The time period chosen has 
no special significance.) The quantity amenable to a stochastic analysis of the Ito-Langevin 



type is the price change during various time intervals [0,|lJ,y,|I^ . To agree with the notation 
used in the previous section we define a; to be a price change during a time interval t. For each 
fixed time interval t we generate a sequence of non-overlapping price changes {xi,X2, ■ ■ .}, 
with Xj = p{jt + 1) — p{{j — l)t + 1). We view the data thus generated for each fixed time 
interval t as a sample selected from a population with some distribution Pmarket{x,t). Thus 
in this application the presumably stochastic variable a; is a price change, and the 'time' 
variable t is really the corresponding time interval. 
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Anomalous diffusion occurs in random systems with correlations in time. Because the 
price change during a longer interval is a sum of price changes during shorter intervals, one 
has 

P{x,ti + t2)= J dXiP{x — Xi,t2\xi,ti)P{xi,ti). (18) 

where P{-\-) is a conditional probability. In the absence of correlation, price changes in 
different intervals are independent, so that P{x — Xi,t2\xi,ti) = P{x — Xi,t2)- In this case 
it is easy to see that the variance grows linearly with time. Anomalous diffusion thus means 
that price changes during successive time intervals are not independent — naturally enough, 
since traders respond to earlier changes. The observed super diffusion of financial markets 
thus indicates correlation, and consequently a nontrivial time dependence. Autocorrelation 
analyses of the S&P500 market have found strong correlations for times on the order of a 
few minutes, with weak persistent correlations at longer times We find the same 



for our data set. Consequently in this work we will concentrate on the intra-day market 
dynamics of intervals less than one hour (1 min < t < 60 min). 

If the dynamics are of Fokker-Planck form, then there must be a consistency between the 
distribution Pmarket{x,t) at a certain time t and the way the distributions evolve in time. 
One would have to find that the data at different times can be fit by distributions of Tsallis 
form, with a Tsallis parameter q independent of time. If so, then the variance will evolve in 
time according to Eq. (§,0) with the same q and appropriate values for D, b. 

We investigate the degree to which the S&P500 price changes at different time intervals 
t can be described by a Tsallis distribution evolving according to a Fokker-Planck equation. 
The approach is to bin the data and perform a direct nonlinear fit of Eq. (^. The 'initial' 
distribution corresponds to the shortest value of t, here to = 1 minute. The = 1 min price 
changes and a fit of Tsallis form are shown in Fig. The fit parameters are q = 1.64 ± 0.02 
and /3{to) = 4.90 ± 0.11, and consequently [from Eq. (D] Z{to) = 1.09 ± 0.02. We then fit 
the data at different time intervals t, with the Tsallis parameter q = 1.64 fixed and f3{t) 
determined by the fit. It is important to this analysis that the data at all times continues 
to be well fit by the Tsallis distribution, and indeed this is what we find. (See Fig. 

The resulting inverse variance (5{t) extracted from the fits is shown in Fig. ^. Time 
evolution controlled by the Fokker-Planck equation predicts the time dependence for (3[t) 
given in Eq. ([T3|). Accordingly we fit this form to the extracted values, finding D = 0.217 ± 
0.003 and b = 0.047 ± 0.004. The Fokker-Planck form for /5(t) is also shown in Fig. (|). 

It is clear from Fig. ^ that the agreement is quite good, indicating time evolution of 
nonlinear Fokker-Planck form. 

As a check on our analysis we also generated simulated data using a Tsallis distribu- 



tion Eq. (^ evolving in time as described in Eqs. (12,15). We then performed the same 
analysis that we performed on the S&P500 data. Briefly, the approach is as follows. We 
use the parameters q,D,b and initial values /3(to), ^'(to) obtained from the S&P500 data. 
Parameters P(t),Z(t) at later times are obtained from Eq. (p!5|,|6|) . At several times t we 
generate a set of random numbers x with the Tsallis distribution Eq. (Rf). This can be done 



by a transformation from a uniformly distributed set of random numbers |T9[, as follows. 
Consider a uniform distribution Puniform(2/) which is unity for < ?/ < 1 and zero otherwise. 
By choosing an appropriate function y{x) we can obtain the desired distribution from the 
uniform distribution in the form 
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PiyX^t^ -Puniform (y) 



dy_ 

dx 



(19) 



This can be solved for y as a hypergeometric function: 



y 



/ P(x',t)dx' = X ■ 2F1 
Jo 



2 g-1'2 



,-,-{q-l)P{t)x'' 



(20) 



For each y generated uniformly, we numerically invert this to find x. The resulting x's are 
distributed according to the Tsallis form, Eq. (^. Thus the simulated price-change data is 
guaranteed to have a Tsallis distribution evolving in time according to the Fokker-Planck 
equation. We then analyzed this simulated data exactly as we did the financial data as a 
check. Fig. p| shows that simulated data closely tracks the actual market data. 



III. CONCLUSIONS 

The nonlinear Fokker-Planck equation for a probability distribution P{x, t) has a time- 
dependent solution equivalent to the distribution obtained by maximization of the Tsallis 
nonextensive entropy. Consequently if a stochastic process has a Tsallis distribution it 
is natural to attempt to work backwards, and ask whether the underlying dynamics are 
Fokker-Planck. We have investigated this possibility in one market, analyzing the high- 
frequency intra-day dynamics of the S&P500 index within this framework. We found that 
the distribution P{x, t) is well fit by the Tsallis form at all times, and evolves in time 
according to the Fokker-Planck equation. Thus the combination of the Tsallis distribution 
and the Fokker-Planck equation unites in a very natural way the power-law tails of the 
distributions and their superdiffusive dynamics. 

It would be of interest to apply this analysis to longer time windows, inter-day (long 
term) dynamics, and also to other markets. 

We acknowledge support from the NSF through grant DMR99-72683. 
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FIGURES 



FIG. 1. Distribution of S&PSOO index price changes at intervals t =1, 10, and 60 minutes, and 
a fit of Tsallis form [Eq. (^)]. (Note: horizontal scales in the three plots are different.) 

FIG. 2. Parameter (3{t) (which is proportional to the inverse of the variance) as a function of 
time interval t. Filled circles: S&P500 index data. Solid curve: Fit of the Fokker-Planck (3{t) from 
Eq. (p^SD- Open diamonds: Simulated data using parameters from the fit. Along most of the curve 
the simulated data are indistinguishable from the actual market data. 
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Fig. 1 F. Michael 
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Fig. 2 F. Michael 



